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ABSTRACT 

We show that there has been a dramatic decline in the abundance of massive galaxies with strong 
stellar absorption lines from 2 : ^ 1.2 to the present. These “HJ-strong”, or HDS, galaxies have 
undergone a recent and rapid break in their star-formation activity. Gombining data from the Gemini 
Deep Deep and the Sloan Digital Sky Surveys to make mass-matched samples (M^ > 10^^-^ Mq, 
with 25 and 50, 255 galaxies, respectively), we find that the fraction of galaxies in an HDS phase has 
decreased from about 50% at z = 1.2 to a few percent today. This decrease in fraction is due to 
an actual decrease in the number density of massive HDS systems by a factor of 2-4, coupled with 
an increase in the number density of massive galaxies by ^ 30 percent. We show that this result 
depends only weakly on the threshold chosen for the HJ equivalent width to define HDS systems 
(if greater than 4A) and corresponds to a (1 + evolution. Spectral synthesis studies of the 

high-redshift population using the Pegase code, treating H^a, EW[0 H], D^dOOO, and rest-frame 
colors, favor models in which the Balmer absorption features in massive H(5-strong systems are the 
echoes of intense episodes of star-formation that faded 1 Gyr prior to the epoch of observation. The 
2 ; = 1.4 — 2 epoch appears to correspond to a time at which massive galaxies are in transition from 
a mode of sustained star formation to a relatively quiescent mode with weak and rare star-formation 
episodes. We argue that the most likely local descendants of the distant massive HDS galaxies are 
passively evolving massive galaxies in the field and small groups. 

Subject headings: Galaxies: formation — galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

The history of cosmic star-formation can be con¬ 
veniently encapsulated by the evolution of the star- 
formation rate density (SERD) as a function of redshift 
2 :. The SERD is a fundamental quantity that places 
models for galaxy formation (e.g. 
In hierarchical models, individual 
galaxies are the ephemeral building blocks of more mas¬ 
sive systems, and galaxy formation is a continuous pro¬ 
cess. In this picture, it is difficult to make an evolu¬ 
tionary connection between a massive galaxy seen at low 
redshift and a very similar galaxy seen at a substantially 
higher redshift, because today’s massive galaxies were in 
several pieces at high redshifts. This difficulty in ‘con- 
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necting the dots’ between similar galaxies over a range 
of redshifts is greatly reduced if massive galaxies form 
earlier and are longer-lived than predicted by early gen¬ 
erations of hierarchical models, as now appears to be 
the case (|Glazebrook et al.ll2004 hereafter Paper HI). 
In this view, a detailed understanding of the formation 
and evolution of individual massive galaxies is entirely 
complementary to the stud y of volume- average d quanti¬ 
ties such as the SERD fe.g. iJnnean et al.ll2005l. hereafter 
Paper V). 

Broad-band colors and spectral features have long 
been used to probe the star-formation history of in¬ 
dividual galaxies, with mixed success. Emission lines 
and blue starlight probe the current star-formation rate. 
Gonstraints on the past history of star-formation from 
the colors and spectral features of older stellar pop¬ 
ulations, however, are hampered by the well-known 
age-metallicity degeneracy. Eortunately, characterizing 
the recent (< 2 Gyr) star-formation history of galax¬ 
ies i s more tractable. In partic ular, the H(5 a Lick in- 
dex llWorthev fc Ottavianm997r and the 4000-A break 
(jDressler fc Shectmanlll987j) . when combined together, 
have proved to be powerful diagnostics for probing recent 
star-formation. Such subtle diagnostics of recent star- 
formation can only be probed by spectroscopy as they 
only have small effects on broad-band colors. Galaxies 
which have experienced a break in their star-formation 
history are particularly well identified by such features. 
Eor instance, the evolution of the H^a Lick index for 
a stellar population reaches high values (> 5 A, cor¬ 
responding to a strong absorption line) only when this 
population is dominated by A stars, which happens ex¬ 
clusively for intermediate ages 300 Myr to ^ 2 Gyr) 


stringent limits on 
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under the condition that no significant star formation 
has occurred in the last few hundred Myrs. This behav¬ 
ior is predicted with good confidence by the recent gen- 
eration of galaxy evolution mode l s, suc h as Pegase.2 


(iFioc fc Rocca-Volmerangd I1997L Il999ll and Pegase- 
HR, (|Le Borgne et al.ll2004ll at low and high spectral res¬ 
olutions, respectively. The systems identified by this par¬ 
ticularly deep H(5 line are usually called “H(5-strong”. 

The 115 line and the 4000-A break have long been used 
to probe stel lar populations, originally mostly in clusters 
of galaxies llDressler fc GunrJ Il98a fCouch fc Sharpleij 


more recentIv also in field galaxies (IZabludoff et al. 

1996 

: IZabludofS 119991: iBaloeh et alJ 

19991: 

Goto et al. 

2003 

: lOiiintero et al.ll2004 iTran et al. 

bioi 

Yang et al. 

2004 

Kaul 

). In a particularly significant recent development, 
Tmann et al.l (1200^1) have used these features to ex- 


amine the star formation histories of large samples of lo¬ 
cal galaxies from the Sloan Digital Sky Survey (SPSS). In 
the pr esent paper we complement the iKauffmann et ahl 
(|2QQ3D work by measuring the same spectral features 
in distant {z = 0.6 — 1.2) galaxies from the G emini 
Deep Deep Survey (GDDS, I Abraham et al.l[2(Kl^ here¬ 
after Paper I). Building two unbiased samples of mas¬ 
sive (M^ > Mq) galaxies at z = 0.6 — 1.2 and 

2 : = 0.05 — 0.1 from the GDDS and the SPSS, respec¬ 
tively, we investigate the differences in the star-formation 
histories at these two epochs. Such a comparison can 
shed light on the evolution of massive galaxies since 
z = 2: the local massive galaxies, dominated by elliptical 
morphologies, are known to be nearly-passively evolv¬ 
ing. Similarly massive galaxies have also been found 
in surprisingly large numbers at higher redshifts (up 
to ^ = 2, Paper III, iFontana et al. |2 004D and some of 
these are already “red and dead” ( Gimatti et al.ll2(K3 
iMcGarthv et al.lT2004 hereafter Paper IV). Our aim is 
to use spectral features to better understand the detailed 
star formation histories of these systems. 

This paper is organized as follows. Our high-redshift 
and low-redshift galaxy samples are described in Sec¬ 
tion [3 Section 13 presents the measurements we have ob¬ 
tained from these data. Results obtained from these mea¬ 
surements are given in Sectional Modeling and inter¬ 
pretation of our data are presented in Sectional and our 
conclusions are given in Sectional Gosmological param¬ 
eters Ho = 75 kms“^Mpc“^, Oa = 0.7, and = 0.3 
are assumed throughout this paper. 


2. SAMPLE SELECTION 

The present study is based on two spectroscopic sur¬ 
veys at low and high redshifts. At < 2 ; >= 1.1, the GDDS 
contains more than 300 high-quality galaxy spectra and 
has a very well-understood selection function (Paper I) 
designed to make it possible to use it as an unbiased 
mass-limited sample. At low redshift {z 2 ^ 0.1), the 
SDSS provides many thousands of galaxy spectra with 
excellent quality. In this Section, we first give a brief 
overview of the common selection criteria for these two 
samples, before describing the details of each sample in 
more detail. 


2.1. Overview 

Our overall sample is selected from the GDDS and 
SDSS using three criteria. (1) We restrict our anal¬ 


ysis to galaxies with secure redshifts and high quality 
spectra. Our threshold corresponds to a typical signal- 
to-noise ratio of S/N > 2 per pixel. (2) Then, we re¬ 
strict the redshift range of each sample to the maximum 
interval over which the measurement of both H(5 and 
the 4000-A break is possible in the data. (3) Finally, 
we select galaxies with a total stellar mass larger than 
IOIO .2 With this purely physical selection we avoid 
color selection biases (for instance, UV-selected samples 
can over-represent star-forming galaxies with respect to 
red galaxies). We choose this particular mass cut-off be¬ 
cause it corresponds to the lowest mass for which the 
GDDS sample can be considered as complete in the red¬ 
shift range defined by our second criterion. Doing so, 
we also maximize the number of galaxies included in the 
high redshift sample, which makes our results more ro¬ 
bust. Finally, given this mass cut-off, we further restrict 
the SDSS sample to the redshift range over which the 
survey is complete. Eventually, we obtain two compara¬ 
ble, mass-limited samples at low (0.05 < 2 ; < 0.1) and 
high (0.6 < 2 : < 1.2) redshifts. 

2 .2. The GDDS sample 

The characteristics of the Gemini Deep Deep Survey 
are described in Paper I; only a very brief description of 
the survey is given here. The GDDS targets galaxies in 
the “redshift desert” (1 < 2 ; < 2). The survey contains 
312 spectra of gal axies in four fields of th e Las Gampanas 
Infrared Survey (|McGarthv et al.|l200lll . The primary 
selection is based on near-infrared photometry {Kg < 
20.6, I < 24.5 in Vega magnitudes) giving preference to 
galaxies with ered I — K color^and photometric^ edshifts 

Spectra were obtained with the GMOS multi-slit spec¬ 
trograph on the Gemini North 8 meter telescope. The 
“Nod & Shuffle” technique was used to optimize sky sub¬ 
traction in order to allow very long exposures (> 30 hours 
per field) to be undertaken. This yielded high-quality 
spectra with FWHM:^ 17 A from which spectroscopic 
redshift determinations were successful for 2 ^ 75% of the 
sample. About 80% of the galaxies with secure redshift 
also satisfy the condition S/N > 2 per pixel around the 
line. 

In this paper our sample is restricted to the 62 galax¬ 
ies with high-quality spectra and secure redshifts in the 
range 0.6 < 2 : < 1.2. This range corresponds to the max¬ 
imum interval over which all galaxies have measurable 
H^a and D^4000 indices. We further limit the sample 
to galaxies with stellar masses > 10^^‘^ M©, leading 
to a subsample of 26 galaxies. The choice of this mass 
cut-off maximizes the number of galaxies in our GDDS 
sample (it is the minimum mass for which our sample 
is mass-complete over the redshift range 0.6 < 2 ; < 1.2, 
as demonstrated in Paper III). Stellar masses, derived 
from template fits to the multicolor (VIK) photometry, 
are taken fro m Paper III. The mas s es use d in this pa¬ 
per assume a iBaldrv fc GlazebrookI llMl^ Initial Mass 
Function (hereafter BG03 IMF). The typical error bar 
on log;L 0 ^^5 estimated by Monte Garlo simulations, is of 
the order of 0.2 dex. A single galaxy (GDDS-02-2197) is 
then removed from this subsample because of very doubt¬ 
ful measured US a and D^4000 indices. This galaxy is an 
outlier with the largest error bars in the sample. In fact, 
these error bars are so large that the object is compatible 
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TABLE 1 

Properties of our sample of GDDS massive galaxies. 


ID 

z 

H5a 

(A) 

EW[0 II] 

(A) 

Dn4000 

logio(M,/M 0 ) g-r SFR([OII]) 
(Moyr-l) 

M,/SFR 

(Gyr) 

weight 

GDDS-22-2863 

0.918 

11.81 ±4.10 

-47.58 ±8.72 

1.16 ± 0.11 

10.30 ±0.30 

0.62 

1.7 

11.6 

0.070 

GDDS-22-0128 

1.024 

11.62 ±2.65 

-56.41 ±12.69 

1.49 ±0.13 

10.55 ±0.28 

0.60 

1.7 

20.7 

0.028 

GDDS-12-8139 

1.189 

8.36 ±0.54 

-6.21 ±0.65 

1.61 ± 0.02 

10.39 ±0.29 

0.40 

2.3 

10.5 

0.019 

GDDS-02-1724 

0.996 

7.58 ±0.19 

-2.45 ±0.32 

1.34 ±0.01 

10.87 ±0.09 

0.37 

2.6 

27.7 

0.007 

GDDS-12-5722 

0.841 

6.14 ±0.58 

-8.74 ±1.11 

1.34 ±0.02 

10.69 ±0.10 

0.67 

1.6 

31.1 

0.166 

GDDS-22-0281 

1.022 

6.10 ±1.26 

-14.88 ±2.31 

1.60 ±0.06 

11.02 ±0.13 

0.71 

1.5 

70.6 

0.166 

GDDS-02-1777 

0.982 

4.74 ±0.90 

-10.81 ±1.71 

1.31 ±0.03 

10.43 ±0.22 

0.56 

2.3 

11.5 

0.033 

GDDS-02-1702 

1.052 

4.58 ±1.58 

-42.99 ±2.42 

1.17 ±0.04 

10.40 ±0.18 

0.33 

9.0 

2.8 

0.020 

GDDS-02-1543 

1.131 

4.22 ±2.65 

-5.34 ±2.71 

1.28 ±0.06 

10.80 ±0.15 

0.67 

0.2 

267.4 

0.070 

GDDS-22-2548 

1.022 

3.47 ±1.51 

-17.08 ±2.46 

1.71 ±0.07 

10.99 ±0.18 

0.70 

2.3 

41.5 

0.166 

GDDS-22-0315 

0.909 

3.13 ±0.96 

-14.90 ±1.05 

1.56 ±0.03 

10.69 ±0.10 

0.53 

4.2 

11.7 

0.013 

GDDS-12-7660 

0.791 

3.04 ±0.79 

-15.66 ±1.11 

1.20 ± 0.02 

10.48 ±0.14 

0.46 

2.5 

12.3 

0.013 

GDDS-02-0715 

1.133 

3.01 ±2.14 

-4.29 ±2.22 

1.35 ±0.06 

11.25 ±0.06 

0.37 

1.8 

101.6 

0.053 

GDDS-15-6851 

1.126 

3.00 ±0.92 

-8.32 ±1.51 

1.64 ±0.03 

10.68 ±0.13 

0.62 

4.2 

11.5 

0.028 

GDDS-12-5337 

0.679 

2.30 ±0.39 

-13.30 ±1.17 

1.67 ±0.02 

10.42 ±0.13 

0.50 

1.9 

14.0 

0.027 

GDDS-12-8983 

0.963 

0.72 ±0.34 

-1.51 ±0.50 

1.77 ±0.02 

10.70 ±0.07 

0.56 

0.6 

89.7 

0.015 

GDDS-22-0893 

0.875 

0.34 ±1.38 

-6.52 ±1.81 

1.93 ±0.07 

10.59 ±0.18 

0.68 

0.6 

60.2 

0.033 

GDDS-02-0857 

1.049 

0.13 ±1.23 

-6.51 ±1.28 

1.88 ±0.05 

11.21 ±0.06 

0.79 

1.8 

87.1 

0.168 

GDDS-12-6800 

0.615 

-0.13 ±0.34 

0.95 ±0.64 

2.03 ±0.02 

10.89 ±0.09 

0.73 

0.1 

1443.0 

0.038 

GDDS-12-6456 

0.612 

-0.43 ±0.57 

-0.06 ±1.08 

1.99 ±0.04 

10.72 ±0.14 

0.73 

0.2 

286.3 

0.027 

GDDS-02-1011 

1.133 

-1.04 ±2.00 

4.93 ±1.53 

1.78 ±0.07 

10.48 ±0.13 

0.65 

0.3 

88.8 

0.070 

GDDS-02-1935 

0.915 

-1.12 ±2.84 

-4.28 ±2.84 

1.82 ± 0.10 

11.12 ±0.07 

0.77 

0.4 

336.7 

0.166 

GDDS-15-5348 

0.964 

-1.97 ±1.19 

-10.82 ±3.55 

1.96 ±0.07 

11.06 ±0.07 

0.71 

0.8 

149.1 

0.166 

GDDS-15-7241 

0.749 

-2.29 ±0.68 

-1.30 ±0.98 

1.74 ±0.03 

10.70 ±0.14 

0.62 

0.2 

260.4 

0.011 

GDDS-22-0510 

0.820 

-2.30 ±1.13 

-2.63 ±1.34 

1.81 ±0.04 

10.51 ±0.13 

0.67 

0.6 

55.0 

0.013 


with the rest of the sample in a formal statistical sense, 
so simply omitting it changes none of our conclusions 
while simplifying the analysis. This leaves us with 25 
galaxies in the final GDDS subsample. The main prop¬ 
erties of these 25 galaxies are given in Table G] A full 
description of all the measured quantities present in this 
table is given in Sectional 

Figure G] presents rest-frame spectra of all 25 GDDS 
galaxies investigated in the present paper, with a 3- 
pixel boxcar smoothing, together with the error spec- 
tra (in blue) and t heir ACS HST images (see also 
I Abraham et al. I l2n?il when available. The color of the 
frame around each image codes the spectral classification 
of the galaxy (Table 5 of Paper I, red for old stellar pop¬ 
ulations, green for intermediate age and blue for young 
stellar populations). The median redshift of this sample 
is z = 0.964. 

The selection function for the GDDS is given in Pa¬ 
per I. The inverse weights in that paper quantify the de¬ 
gree to which each galaxy is over-represented in di I — K 
vs. / color-magnitude diagram. It is important to note 
that the primary goal of the GDDS survey was to study 
red galaxies at high redshift. Therefore, by design, this 
survey specifically targeted red galaxies with a sparse 
sampling rate of ^1/2 compared to ^1/7 for the remain¬ 
der of the photometric sample. The inverse weights must 
therefore be used to account f or th is unbalance. As we 
have just described in Section [TTl in the present paper 
further restrictions on the sample are introduced (so as to 
only include galaxies with high confidence redshift class 
and high spectral quality), and a new set of weights is 
needed to define the selection function for this sample. 


These are given in the last column Table . In Sec¬ 
tion ^ we will show that our results are not strongly 
dependent on the adopted weights. 

2.3. The SDSS sample (z 0,1) 

At low redshift, we define a sample o f galaxies drawn 
from the Sloan Digital Sky Survey DR2 (|Abazaiian et al.l 
1200411 . We restrict the present investigation to the ‘main’ 
galaxy sample, i.e. galaxies with Petrosian r-band mag¬ 
nitudes in the range 14.5 < r < 17.77. Only galaxies 
with secure redshift measurements are considered. To re¬ 
duce sampling incompleteness, we focus on galaxies with 
2 : > 0.05, and in order to be consistent with the quality 
of the GDDS spectra, we also exclude the few galaxies 
(0.2% of the sample) with S/N < 2. 

As for the GDDS sample, we also impose a mass 
cut, so that our sample includes only massive galaxies 
(M^ > 10 ^^-^ Mq). Ac cording to the SPSS stellar mass 
function determined bv iPanter et al.l (|2004ll . the sample 
is complete for \ogiQ{M^/h~‘^) > 10.2 only at 2 ; < 0.1. 
Using the cosmology adopted in our study, this corre¬ 
sponds to Howeve r, this mass-function 

is estimated using a lSaloet^ (|1955fl IMF, while the stel¬ 
lar masses of GDDS galaxies are estimated using the 
BG03 IMF. The conversion factor between these mass 
estimates is Msp = 1.82Mbg (Paper III), or 0.26 dex. In¬ 
cluding this correction, we consider that the SDSS sam¬ 
ple is complete for Mq (using the BG03 

IMF) for 2 ; < 0 . 1. W e use the mass estimates from 
iKanffmann et al.l (|2003fl which are corre cted from dnst- 
obscuration. Because these authors use a lKronoal ()200lD 

We emphasize that these are inverse weights. A small number 
implies that the object needs to be weighted more heavily because 
it is under-represented in the sample. 


















4 


D. Le Borgne et al. 


GDD5-ZZ-Z8b3 ■ ... . 

z=o.9i8 I :: : : 


[oil] KH H6 Hy 


GDDS-22-0128 
z= 1.024 


.hit. 



ODDS-15-6851 

z=1.126 

H6*=3.00 


ODDS-12-8139 

z=1.189 

H6a=8.36 



ODDS-12-5337 

z=0.679 

Ha*=2.30 


GDDS-02-1724 

z=0.996 

H6*=7.58 





010 

GDDS-12-8983 
z=0.963 

H6a=0.72 




001 

GDDS-12-5722 
z=0.841 

H6*=6.14 

V —-^- 





110 

/ 

GDDS-22-0893 

z=0.875 

H6*=0.34 




001 

hrWWVY 

—TvAa— 

GDDS-22-028i 

z= 1.022 

H6*=6.10 




*Mv/ 

. '"oil 

. 

GDDS-02-0857 
z= 1.049 

H6*=0.13 


1 


A 1 001 


o< 


O 

CO 

\ 

cn 

0 

X 

D 


TD 

0) 


GDDS-02-1777 
z=0.982 

H6*=4.74 j 

IfJV^ 




010 

# 

GDDS-12-6800 
z=0.615 

H6*=-0.13 




001 

GDDS-02-1702 1 

z= 1.052 

H6^=4.58 




"—\ 

110 

0 * 

GDDS-12-6456 
z=0.612 

H6^=-0.43 





O 

o 

GDDS-02-1543 

z=1.131 

H6^=4.22 


fv 

Iffy 


001 

GDDS-02-10li 

z=1.133 

H6^= -1.04 



lU'V' 


001 

4» 

GDDS-22-2548 
z= 1.022 

H6*=3.47 1 





GDDS-02-1935 

z=0.915 

H6a=-1.12 




001 


GDDS-22-0315 

z=0.909 : : 

H6a=3.13 J 




011 

• 

GDDS-15-5348 
z=0.964 

H6a=-1.97 




vfl 

001 

GDDS-12-7660 
z=0.791 

H6a=3.04 I 





110 

0 

GDDS-15-7241 
z = 0.749 

H6*=-2.29 





■ 001 

‘ t 

GDDS-02-0715 ' ' 

z=i.i33 : : 


/iljlfW 

k 

110 

% 

GDDS-22-05i6 

z=0.820 

H6a=-2.30 





o 

o 


5000 3500 4000 ^ 4500 

Wavelength (A) 


3000 3500 4000 ^ 4500 

Wavelength (A) 


Fig. 1.— Sample of the 25 mass selected GDDS galaxies, sorted by decreasing H(5 a (expressed in A). The error spectra are in blue and 
HST ACS images are shown when available. The size of these images is 5x5 arcsec^. The color of the frame around each image codes the 
luminosity-mean age of the galaxy (estimated by visual inspection of the spectrum and broad-band SED, cf Paper I): red for old stellar 
populations (codes 001 and Oil), green for intermediate age (code 101) and blue for young stellar populations (code 110). Note that slightly 
negative values for H(5 a (;$ SA) do not necessarily correspond to emission lines (e.g. GDDS-15-7241), because of the very definition of this 
Lick index: even when the H(5 line is in absorption, the pseudo-continuum can be “lower” than the line itself. 


IMF, we subtract 0.07 dex from their masses to convert to 
BG03 IMF mass estimates. For our subsample of SDSS 
galaxies, the typical 95% confidence range on logM^ is 
0.2 dex. 

After imposing the cuts described above, our final 
SDSS low-redshift sample contains 50, 255 galaxies. 


In this Section, we provide a detailed description of 
the quantities measured from our low-redshift and high- 
redshift samples. The reader who is not interested in de¬ 
tails of how these measurements were made may wish to 
simply note that we measured a single rest-frame broad¬ 
band color and the spectral indices listed in Tabled and 
then move ahead to Sectional 


3. MEASUREMENTS 
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TABLE 2 

Definitions of the spectral indices used in this paper. 



Blue band 

Feature band 

Red band 

Index 

(A) 

(A) 

(A) 

D„4000 

3850.00-3950.00 


4000.00-4100.00 

H5a 

4042.85-4081.00 

4084.75-4123.50 

4129.75-4162.25 

EW(H(5) 

4017.00-4057.00 

4083.50-4122.25 

4153.00-4193.00 

EW[0 II] 3696.30-3716.30 

3716.30-3738.30 

3738.30-3758.30 


3.1. Rest-frame g-r eolor 

At low redshift, SDSS galaxy spectra and fiber u,g,r,i,z 
magnitudes are publicly available. However, given that 
we use the SDSS spectra to measure other spectral in¬ 
dices (D^dOGO and we decided to measure a color 

directly on the SDSS spectra to be consistent in our 
measurements (and to avoid, for example, aperture ef¬ 
fects which might affect differently the spectral features 
and the photometry). Because of this choice, g — r and 
B — V are the two main optical colors available for the 
SDSS galaxies. We checked that both colors give simi¬ 
lar results, and we concentrate on the g — r color in the 
following. The g — r color is measurable on 95% of the 
SDSS spectra. 

Unfortunately, because of the limited wavelength cov¬ 
erage of the GDDS galaxy spectra, it is not possible to 
measure a standard optical color from them in the same 
way. At redshift 1 . 1 , the optical spectra of the GDDS 
galaxies correspond to UV rest-frame spectra. However, 
the existence of optical and near-infrared photometry 
makes it possible to estimate a rest-frame optical color, 
using k-corrections, in the following way. The rest-frame 
g — r color of the GDDS galaxies was measured on the 
synthetic spectra that fit best the observed spectral en¬ 
ergy distributions, using the B,V,R,I,z,J,H and K pho¬ 
tometry. The procedure used to fit the data combines 
two instantaneous bursts of star formation at different 
ages and metallicities, reddened by dust (parameterized 
by E(B-V), ranging from 0 to 2 ). A minimization of the 
using the ages, the metallicities and the extinction as 
free parameters, provides a best-fitting synthetic spec¬ 
trum. (It is worth mentioning that this procedure was 
not used to derive ages, metallicities, or extinction pa¬ 
rameters: doing so would require a much more refined 
grid than the one we used because of the degeneracies 
that exist in the measurement of these three parame¬ 
ters.) Column 8 of Tabled lists the measurements of the 
g — r color for the GDDS galaxies in our sample. 

3.2. D^4000 

The g — r color of a galaxy can be affected by the pres¬ 
ence of dust, and the D^4000 at 4000 A is often used 
to estimate ages instead because it is much less sensitive 
to extinction. For example, a reddening of E(B-V)=0.5 
shifts D^4000 by less than 10 percent, which is of the or¬ 
der of the change produced by a relative difference in age 
of 0.25 dex (or 45%). This effect is small in compar¬ 
ison to the corresponding reddening of the g — r color, 
which can be mistaken for a 1 dex age effect. Re¬ 
markably, the dependence of D^dOOO on the metallicity 
is much stronger than its dependence on the extinction: 
for a 10 Gyr old population, it varies from 1.7 to 2.7 
when the metallicity varies from Z^/b to 2 . 5^0 (e.g. 


Kanffmann et al.ll2003fl . The definition of D^4000 (from 
Balogh et al.lll999ll is reported in Table [2| together with 
the other indices used in this paper. 

3.3. The H5 line 

In the absence of emission, the equivalent widths of 
Balmer absorption lines decrease steadily with time for 
an old (age> 1 Gyr) passively evolving stellar popula¬ 
tion. Moreover, they are quite generally little to the 
mean stellar metallicity and dust content of a galaxy. 
However, in practice, Ha and H/^ are often partially filled 
in by nebular emission, making their interpretation more 
difficult. Very high order lines are generally too weak to 
be easily measurable, so H 7 and Hii are the two lines 
of choice for probing stellar ages. In the present study 
we concentrate on 116 because, unlike the lower order 
Balmer lines, it is measurable in optical spectra of galax¬ 
ies redshifted up to ^ = 1 . 2 . 

An interesting property of the BS absorption line lies 
in its behavior between the birth of a generation of stars 
and the age of 2 Gyr: the equivalent width first increases 
strongly during a few 10^ years, reaches its maximum 
value (almost 10 A) when the population is dominated 
by A stars, and then decreases continuously, passing 
through its initial value when the population is about 
2 Gyr old. This behavior implies that high equivalent 
widths in absorption (> 4 A) can only be reached in 
galaxies where some substantial star formation occurred 
a few 10^ years before, and where the subsequent star 
formation events were weak. We will refer to systems 
with 116a > 4 a as being “H(5-strong” (HDS). 

To measure the equivalent wid th of HJ, we use the 
defin ition of the Lick index HA llWorthev fc Ottavianil 
I 1997 II which is well suited to probe stellar populations 
dominated by A stars in the optical. We degrade the 
resolution of all the rest-frame spectra to the appropri- 
ate Lick resolution (1 1 A), following the prescription of 
iWorthev fc Ottaviani lll997li. 

For the sake of consistency between the two data sets, 
we choose to measure the BSa and D^4000 indices our¬ 
selves on the SDSS spectra, rather than using the pre¬ 
viously published values. To validate these measure- 
ments, we compare the m with the values published by 
iKanffmann et al.l (|2003L hereafter K03) in Fig. (21 Perfect 
agreement is not expected: the differences between our 
measurements of BSa and the values given by K03 can 
be partly explained by the different spectral resolutions 
(:^ 3 A for K03, and ^ 11 A for this work which uses 
the Lick resolution). Moreover, unlike those authors, we 
choose not to correct BSa from the nebular emission for 
two reasons: first, we prefer to compare the raw data to 
models that include the contribution from emission lines; 
second, the correction from the contribution from emis¬ 
sion lines to the H^a index is too uncertain for the GDDS 
galaxies which have on average lower S/N ratios than the 
SDSS spectra. However, the agreement between the two 
measurements is good, given the different considerations 
involved in the measurements. For large values of Hi^a, 
our measured values are smaller than the K03 values by 
~ 1 A. This is to be expected as a correction for line 
emission in such young populations would increase BSa 
by about 1 — 3 A, as predicted by the Pegase models . 
A similar correction is mentioned inlBalogh et all (|l999ll , 
and is confirmed bv iGoto et all (|2003fl who measured a 
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Dn4000 (K03) (K03) 

Fig. 2. — Compariso n of Drt4000 and H(5a m easurements on our 
SDSS sample between IKauft'mann et all 1200,^111 (with a correction 
for nebular emission and a higher spectral resolution) and this work 
(including the emission lines and at the Lick spectral resolution). 
The grey scale level represent the number of galaxies in each pixel 
in a logarithmic scale, darker pixels having more galaxies. 


median correction of 15% for the 116 equivalent width on 
a sample of galaxies drawn from the SDSS. 

3.4. Emission lines 

The nebular emission lines produced in H II regions 
are excellent indicators of recent star-formation. In the 
redshift range of the GDDS sample, the main visible lines 
are [OIIJA3727 and the 116 Balmer line (Hq and Up are 
too much redshifted to be visible). In the following, we 
concentrate on the equivalent width of the [OIIJA3727 
line as an indicator of ongoing or very recent (less than 
50 Myr) star formation. 

3.5. Errors 

The uncertainties on the measurements of H(5 a and 
D^4000, and E W[0 II], repor t ed in Table ^ are com- 
puted following iCardiel et all 1(12^, using, for each 
galaxy, the measured mean signal-to-noise ratios inside 
the integration bands of Table [2j To confirm these er¬ 
ror bars, Monte-Carlo simulations were carried out as¬ 
suming a Gaussian distribution for the error spectrum. 
For each galaxy, 50 simulations were considered, and the 
rms dispersion of the distribution of indices matches very 
well the erro r bars predicted using the formulae from 
ICardiel et'all lll998ll : the rms for the distribution of rel¬ 
ative differences between these two errors is 0.22 for 116a 
and 0.12 for D^4000. 

Although we believe that the observed fractions of HDS 
galaxies in each redshift range are robust, several poten¬ 
tial sources of errors should first be considered. 

The very definitions of the indices used might be a 
source of systematic errors. We confirmed that the use 
of the EW(H(5) definition instead of the index 116a does 
not change any of our conclusions. 

Some errors may also come from the calibration of 
the original data. In particular, the reddest end of the 
GDDS spectra is known to be subject to a quite un¬ 
certain correction (the ”redfix” correction, cf. Paper I) 
due to charge diffusion in the GCD detectors in the red. 
Fortunately, this correction applies to the portion of the 
rest-frame spectrum which contains the 116 line only at 
redshifts larger than 1.3. Therefore, it does not affect our 
study. The selection function, which introduces impor¬ 
tant weighting, might also be a concern. However, our 
results are qualitatively the same if this selection function 


is omitted. Finally, we estimate that the aperture effects, 
which could potentially affect the observed integrated 
spectral indices and colors, are small on average. Indeed, 
as GDDS galaxies are observed at higher and higher red¬ 
shift, a larger fraction of the galaxies in enclosed in the 
0.75” aperture slit (corresponding to a physical diameter 
of 4.7 kpc at 2 : = 0.6 and 5.8 kpc at 2 : = 1.2), and more 
light from the outer parts contributes to the integrated 
indices. However, over the whole 0.6 < 2 ; < 1.2 inter¬ 
val, a larger fraction of the total light is measured than 
that measured by the SDSS fi bers on the ^ = 0»Q 5 — 0.1 
gal axies. The recent wo rk of iBrinchmann et al.l (120041) 
an d iKewlev et al.l (|2004D showed that aperture effects for 
emission-line-derived SFR in the SDSS galaxies is negli¬ 
gible for 2 : > 0.05. The aperture effects on the stellar 
absorption lines (typically more concentrated than the 
emission-lines) are expected to be negligible. 

Overall, the main sources of uncertainties are unlikely 
to explain the observed trend of the decrease of the frac¬ 
tion of HDS galaxies since 2 ; = 1.2 described in next Sec¬ 
tion. We conclude that this evolution must be real and 
can only be explained by changes in the star formation 
histories of massive galaxies during the past 10 Gyr. 

4. COMPARED PROPERTIES OF LOCAL AND DISTANT 
MASSIVE CALAXIES 

4.1. General results 

Measurements of the colors and spectral features in the 
galaxies in our GDDS sample are presented in Tabled 
Figures ini and 0 show the distributions of H(5 a, D^dOOO, 
and g — r rest-frame color as a function of redshift. In 
these figures, GDDS galaxies categorized as H(5-strong 
(H^a^ 4 A) are denoted with open circles, and the re¬ 
maining galaxies of the GDDS sample are shown as filled 
circles. In addition, these figures show two-dimensionally 
binned distributions (for both the GDDS sample and the 
SDSS sample) using colored cells. The reddest (or dark¬ 
est) cells contain the largest number of galaxies. The 
size of the cells is chosen to be on the order of, or larger 
than, the typical error bar of the various observed quan¬ 
tities. The density of GDDS galaxies in each binned cell 
has been corrected using the weights given in Table Q 
and Monte-Garlo simulations of each GDDS galaxy were 
used to estimate the true distributions shown in these 
figures. (Note that this explains why the spread data in 
the color-coded cells is a little larger than the coverage of 
the individual points). Each of these figures also shows 
evolutionary models which will be described in the next 
section. 

A comparison of the distributions of the GDDS and 
SDSS samples in these figures shows a number of inter¬ 
esting results. The most striking difference between the 
low-redshift and the high-redshift samples is the dispro¬ 
portionately large number of H(5-strong systems at 2 ; ^ 1 
in the GDDS sample, clearly visible in Figure 01 

To quantify this observation, we measure the number 
space densities (per comoving Mpc^) of massive galaxies 
and massive HDS galaxies in various redshift intervals. 
For GDDS galaxies (0.8 < 2 ; < 1 and 1 < 2 ; < 1.2), we use 
the sam pling weights d efined above and the 1/Vmax for¬ 
malism (|Schmidtil 19681) to correct from the fact that the 
faintest galaxies in our sample could not be observed if 
they were at higher redshift. At low reds hift, we use the 
SDSS sample and the Vmax values from iBlanton et alJ 
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Lookbock time (Gyr) 

0123456 7 8 9 



Fig. 3.— Evolution of the H(5a index from 2 ; = 1.2 (GDDS sample) to the present (SDSS sample). For each sample, the plane is gridded 
and the number of galaxies that can be counted in each cell (we call this quantity “density”) is color-coded: a density of 1 (reddest or darkest 
color) is assigned to the cell containing the most objects. The GDDS densities are converted from raw number counts by using the selection 
function and Monte Garlo simulations. The GDDS individual galaxies are represented by filled (H(5a< 4 A) and empty (H(5a> 4 A) circles. 
The models described in Section |3 are over-plotted: “burst” mode models are represented by black lines and “quiescent” mode models by 
red lines. The two ’cycle’ models plotted in the following figures are omitted here for clarity: they would blur the figure, because their H(5 a 
indices alternatively reach very positive (when star-formation has ceased) and very negative values (because of emission lines) on short 
time-scales. 


(|200dl) . At every redshift, we treat separately the pop¬ 
ulation of HDS galaxies, defined by 4 A. Fig¬ 

ure 0 shows the evolution of the number densities and 
mass densities from z ~ 1 to ^ ~ 0 for massive galax¬ 
ies (M^ > Mq with the BG03 IMF). In par¬ 

ticular, it shows that the z ^ 1 massive galaxies are 
about one third (in number) of the local massive galax¬ 
ies. This means that the GDDS galaxies in our sam¬ 
ple are almost certainly progenitors of some of the local 
massive galaxies (because the stellar mass is expected 
to either increase or keep steady), but about two thirds 
of the local massive galaxies had progenitors at z ~ 1 
that were not massive enough to pass the mass thresh¬ 
old of Mr:^. The figure also confirms prev i- 

ous results (|Glazebrook et al1l2004 iFontana et al.ll2004D , 
showing that the total stellar mass contained in mas¬ 


sive galaxies at in ^ 1 is already about a third of its 

present-day value . It is more diffic ult to link this result 
with the claim of iBell et al.l (|2004D that the stellar mass 
in red galaxies must have doubled since z = 1. Indeed, 
our mass-selected sample contains not only galaxies on 
the red-sequence, but also star-forming galaxies (see Ta- 
bleP). Moreover, the mass cut-off adopted here probably 
excludes the faint-end of the red-sequence dit z ^ 1. 

It also appears for the first time that although the 
number and mass densities of all massive galaxies both 
increase by a factor of 2-3 from 2 ; = 1 to 2 : = 0, these den¬ 
sities follow the opposite trend for HDS massive systems: 
they actually decrease by a factor of 2-4 (for number den¬ 
sities) and 2-7 (for mass densities) during the same pe¬ 
riod. As a result, the fraction of massive HDS galaxies is 
indeed much higher dit z ^ 1 1/2) than a.! z ^ 0 (less 
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Fig. 4. — Similar to Fig. I^but for two other tracers of the luminosity-weighted age of galaxies: Dn4000 and the g — r rest-frame color. 
Symbols as in Fig.|3 


than 1/10). A direct consequence of this strong decrease 
of number densities for HDS galaxies is that most of the 
z = 1 HDS massive galaxies necessarily end up in lo¬ 
cal non-HDS (and probably passively evolving) galaxies. 
This disconnect between distant and local HDS galaxies 
is not too surprising because the timescales involved are 
quite different: less than 1-2 Gyr for the HDS phase and 
^ 7 Gyr for the Hubble time interval between these two 
redshifts. 

In the second part of this section we will consider how 
best to parameterize this evolution, and explore the ex¬ 
tent to which the evolution depends on the equivalent 
width threshold beyond which galaxies become flagged 
as HDS systems. 

The comparison of SDSS and GDDS samples in Fig¬ 
ure 0 shows complementary information. The non-HDS 
GDDS galaxies (black dots) tend to have D^dOOO indices 
almost as large as the values observed at low redshift, 
whereas the HDS galaxies (empty circles) have a smaller 
D^4000 index betraying a smaller luminosity-weighted 
age. This distinction between HDS and non-HDS galax¬ 
ies is much less clear for the g — r rest-frame color (right 
panel of Fig. which can be attributed to differences 
in extinction by dust. However, over the redshift interval 
explored by this paper, the GDDS sample is bluer overall 
than the SDSS sample. 

Figure |S| shows, for each sample, the correlation be¬ 
tween H(5 a and D^dOOO. As will be described in the 
next section, the distributions shown in this figure clearly 
demonstrate that both the SDSS and GDDS samples 
of massive galaxies can be divided into two categories: 
old evolved galaxies with a large D^4000 break and a 
small H^a index, and galaxies dominated by younger 
stellar populations, with a smaller D^4000 break betray¬ 
ing some recent star formation activity. The well-known 


bimodality in the pro perties of blue and r ed populations 
of SDSS galaxies (e.g. iBaldrv et al.ll20d3l is only faintly 
visible in the left panel of this figure; bimodality is not 
particularly striking in our data because our mass selec¬ 
tion isolates the reddest part of the bimodal distribution. 
Both SDSS and GDDS samples of massive galaxies con¬ 
tain a large fraction of old, passively evolved galaxies 
(the reddest cells in Figure EJ. It is once again seen 
that, amongst galaxies that have recently formed stars 
(Dn4000< 1.5), the fraction of HDS galaxies is larger in 
the GDDS than in the SDSS. Another interesting infer¬ 
ence regarding the physical properties of the HDS mas¬ 
sive galaxies can be drawn from Figure El most of the 
HDS galaxies seem dominated by young stellar popula¬ 
tions (they all have a small D^4000 break). 

Figure 13 presents correlations between }ISa and g — r 
rest-frame color. Most of the HDS objects have redder 
g — r colors than predicted from passively evolving mod¬ 
els (which include no or very little dust). This suggests 
that massive HDS galaxies, from 2 : = 1 to the present, 
share a common property: they are considerably more 
reddened by dust (E(B-V)=0.2 —0.5) than evolved galax¬ 
ies. This inference is perhaps not too surprising: hav¬ 
ing experienced recent episodes of intense star formation, 
they are more likely to be dusty. This property is con- 
fir mnd__by_^MwoMm]ni^oi^^ test, implemented as 

in iFasano fc Franceschinil (|l987f) . The two-dimensional 
statistic D is 0.9, and a Monte-Carlo simulation ensures, 
at a > 99.9% confidence level, that the data is not con¬ 
sistent with dust-free models of HD S galaxies. This con- 
clusion is in good agreement with iShiova et al.l (|2004ll 
who suggest an Ay >0.5 mag extinction in the red HDS 
systems. It is worth noting that the effect of dust on the 
spectral indices is very weak: our tests show that E(B- 
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^ GDDS (this work) 

Blanton et ol. 2003 
Q Cole et ol. 2001 
^ Glozebrook et ol. 2004 . 



1.0 

Redshift 


2.0 


Fig. 5. — Evolution of the number and mass densities of massive 
galaxies and massive HDS (H(5a> 4A) galaxies. All these measure¬ 
ments concern only galaxies more massive than Mq. 

Sampling weights are applied to GDDS galaxies, and the l/Vmax 
formalism is used at every redshift. Poisson error bars are shown 
for the z ^ 1 GDDS galaxies. These are too small to be shown 
for the z 0 s ample s. Densities at low redshift were taken from 
IBlgmton et al.l I200ill) (from SDSS data, crosses) and IGole et ahl 
11200111 (from 2dF data, squares), and corrected for the same IMF 
and c osmology. Mass densities from Paper III IGlazebrook et ahl 
l2(^ are shown in the lower panel for reference with a black solid 
line; the corresponding number densities are shown in the upper 
panel. The high redshift (GDDS) and low redshift (SDSS) samples 
used in this study are linked with dashed lines. 


V)=0.5 would not increase 116 a by more than 0.2 A and 
D^4000 by more than 10 percent. 

4.2. Evolution of the fraction of HDS systems 

Figure [S| presented evidence for a differential evolution, 
from 2 ; = 1.2 to the present, of the number and stellar 
mass densities of HDS massive galaxies with respect to 
populations of massive galaxies. This evolution can also 
be measured and parameterized more simply through the 
fraction /hds of HDS galaxies in unbiased samples of 
massive galaxies at various redshifts. 

Figure |HI presents the evolution of this fraction and 
clearly shows a steady increase of /hds with look-back 
time, going from /hds= 0-01 ~ 0-2 in the local universe 
to /hds— 0-4 — 0.8 at 2 : = 1.2. The different panels 
in this figure illustrate cases where the H^a threshold 
(used to define HDS systems) is H^Amin=3,4, and 5.5 A, 
and show that the measured fraction of HDS systems in¬ 
creases markedly in all cases. In this section, we propose 
ourselves to parameterize this evolution as a function of 
the look-back time and as a function of the threshold 

H^A min- 

But first, it is worth investigating how sensitive this ba¬ 
sic result is to the the underlying sample selection func¬ 


tions of the SDSS and GDDS. To investigate this, we ex¬ 
plore the effect of the different selection functions applied 
to the datasets. For the SDSS sample, which is effectively 
complete, we can simply measure the fraction of HDS 
galaxies in a single small redshift bin: 0.05 < 2 : < 0.1. 
For the GDDS sample, we define three redshift bins with 
Az = 0.2 between 2 : = 0.6 and 2 : = 1.2 and we mea¬ 
sure the fraction of HDS galaxies in each bin, weighting 
the contribution from each individual galaxy by the sam¬ 
pling weights listed in Tabled In effect, each galaxy is 
treated as if it were “1/weight” galaxies. The effect of 
applying these weights on the measured /hds 5 and of ig¬ 
noring them completely, can be assessed by comparing 
the diamonds and circles in Figure |H1 Glearly the effect 
of the weighting by the selection function is modest. This 
is probably due to the fact that the vast majority (23/25) 
of the massive galaxies in our sample have an observed 
I-K color greater than 3 magnitudes and lie in a region 
of the I vs I-K diagram where the GDDS sampling effi¬ 
ciency is fairly high (cf Fig. 12 and 13 of Paper I: about 
half of the red galaxies were targeted and a redshift could 
be measured for most of them). 

It is worth noting that the error bars on the points 
shown in Figure |HI cannot be properly described by strict 
Poissonian statistics, because the quantities plotted are 
ratios. Instea d, we compute the uncertainties following 
the method of lPatern^ (|2004f) which is based on the ap¬ 
plication of the Bayes’ Theorem to binomial statistics. 
This method gives realistic errors even when the total 
number of galaxies inside a bin is very small (and the 
error on the ratios computed via Poisson errors break 
down), and when the fraction is close to 0 or 1 (where 
strict binomial errors break down). We choose to repre¬ 
sent the shortest 68.3% confidence interval^^. Another 
precaution is necessary: because the error bars for the 
H^a index are quite large for the GDDS sample, some 
galaxies near the limit adopted for H^a could be consid¬ 
ered as HDS or not. To quantify this effect, we again 
made a Monte-Garlo s imul ation for each galaxy, exactly 
as described in Section lTTl For each simulation, we mea¬ 
sure the fraction /hds and the associated lower and up¬ 
per confidence limits. The fraction that we finally adopt 
is the median value of the distribution of /hds 5 and the fi¬ 
nal lower and upper limits are taken as the 68 % quantiles 
of their respective distributions. The final confidence in¬ 
terval is of the same order of the confidence interval when 
only binomial statistics are considered. 

In addition to showing clear evidence for evolution. 
Figure IHI also shows that the choice of the threshold for 
H^a does not strongly influence the basic trends seen. 
We will now show that it is quite straightforward to pa¬ 
rameterize /hds as a function of both redshift and of the 
threshold H^Amin- 

At a fixed redshift, /hds obviously depends on the 
threshold H^Amin- The top panel of Figure IHI presents 
the histogram of the distribution of HSa for the SDSS 
galaxies. It is very well fitted by an exponential dis¬ 
tribution (a straight line in this log-linear diagram) for 
H(5 a^ 4 A: at the median redshift of 2 : = 0.075, we have 

logio [clN(z)/dmA] = A(z) + a(z) ^A (1) 

The confidence limits derived using this method are compat¬ 
ible within 4 % with the lim its computed from the approximate 
expressions of lGehrel^ 1198611 . 
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Fig. 6. — Variation of H(5a versus 0^4000 for SDSS galaxies (left panel) and GDDS galaxies (right panel). Symbols as in Fig. |3 The 
arrows on the models (lines) indicate the direction of the evolution. For the ’cycle’ models (closed loops), only the steady-state cycle (at 
2 ; = 0) is shown, for clarity. 



g-r 


g-r 


Fig. 7.— Variation of H(5a versus g — r rest-frame color for the 
as in Fig. 1^ 


with a = —0.52. 

We assume in the following that the form of this dis¬ 
tribution is still valid at higher redshift. We will show 
that the data available at 1 is compatible with this 
assumption. 

A distribution of this form should naturally lead to a 
linear equation for log^g [/HDs(H(5Amin)] at a fixed red- 
shift: 

logio [/HDs(2,H(5Amin)] = Oi {z) [H(5a min “ H(5 a(2:)] 

n (2) 

with a' = am the ideal case. is introduced here 


sample (left panel) and the GDDS sample (right panel). Symbols 


to account for the offset of this linear equation at the 
redshift 2 :. 

The top panel of Figure 11171 shows this quantity /hds ia 
the SDSS, as a function of the threshold Hi^Amin- As ex¬ 
pected, /hds changes with Hi^Amin following equation (21 
with a' = —0.53 ~ a. 

At higher redshift, if the distribution of 116a is still a 
exponential one, then the fraction of HDS galaxies should 
also follow equation [3 The three bottom panels of Fig¬ 
ure 13 present the histograms of the distributions of H^a 
for GDDS galaxies in the three redshift bins with mean 


SDSS 
































































GDDS VI - Massive H(5-Strong galaxies at ^ 1 


11 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
Redshift Redshift Redshift 


Fig. 8. — Evolution with redshift of the fraction /hds of H(5-strong galaxies in populations of massive galaxies. Each panel corresponds 
to a different threshold for the minimum H(5 a value defining HDS galaxies. Triangles represent the SDSS sample. Diamonds and circles 
represent measurements on the GDDS sample with and without weights, respectively. The increase of /hds with redshift (roughly, 
logio(/HDs) oc — (1 + corresponds to a (1 + z)^-^ evolution in the distribution of the H(5 a index. The dotted lines illustrate this 

evolution parameterized by equations |2]and|3 (see text for details). 


values z = 0.7,0.9 and 1.1. These distributions are the 
results of the Monte-Carlo simulations presented in Sec¬ 
tion^ and have an arbitrary normalization. As for the 
low redshift sample, they seem to be reasonably well de¬ 
scribed by exponential distributions, with the remark¬ 
able point that the slope a{z) increases with redshift. 
A graphical measure of this increase of a{z) is given in 
Figure ITTl (top panel). We choose to parameterize this 
evolution by an expression of the form 

a{z) = k{l^z)-^ (3) 

and we find that m = 2.0 ±1.4 and k = —0.6. 

As we said before, if the distribution of H^a remains 
close to a exponential one at high redshift, the direct 
measure of /hds as a function of 116 a min should evolve 
following equation 121 In other words we should find that 
a'(z) ~ a{z). The two lower panels of Figure ITOI present 
our measurements of /hds as a function of 116 a min for the 
two highest redshift bins. We observe that at the three 
redshifts z:^0.1, 0.9 and 1.1, /hds indeed follows Equa¬ 
tion (21 with ~ 1.5 A in every case. (An explanation 
for the interesting property that this value seems to be 
a constant is reserved to further investigation.) Com¬ 
paring with the previous figure, we do indeed see that 
a'(z) ^ a(z) in these three redshift bins. FigurelTTKbot- 
tom panel) shows the evolution of o:'{z), and the best fit 
is obtained for m' = 2.5 ± 0.7 and k' = —0.65, if we use 
the same form as in equation 121 The values of m and 
m' are compatible within their error bars but they are 
slightly different. We can think of two reasons for this 
(small) difference. Either the distributions of USa are 
not really exponentials at high redshift, or the difference 
is simply due to the smallness of our sample. A larger 
sample would be useful to refine these values. In any 
case, if the distribution is indeed a exponential, then we 
believe that m', derived from the fits of /hds(^, Hc^Amin), 
is actually a more robust number to measure than m, 
simply because it measures a cumulative quantity. The 
best-fit evolution of /hds (with the parameters m' = 2.5, 
k' = —0.65 and H^a = shown explicitly in Eig- 

ure IHI (dotted lines) and matches the data quite well. 


One might notice that the data point for the redshift bin 
2: = 0.6 — 0.8 is slightly below the best-fitting idealized 
curves: in this particular redshift bin, our sample con¬ 
tains only four galaxies, and this small number might not 
be statistically very significant. Moreover, the proposed 
parameterization should only be valid for Hi^Amin^ 4 A, 
and the left panel of Eigure|H|is outside this range, with 
Il(5Amin=3 A. Therefore, a good agreement is not ex¬ 
pected in this case. 

To summarize, we have shown that the fraction of HDS 
galaxies can be parameterized by equations (21 and 01 
/hds depends both on redshift and on the threshold 
H^Amin, but thcsc two Variables are independent. In par¬ 
ticular, we stress that the evolution of /hds with redshift 
does not depend on the threshold 116a min chosen to de¬ 
fine HDS galaxies, as long as Hi^Amin^ 4 A. We have also 
shown that this evolution of /hds is a direct consequence 
of an evolution of the RSa distribution. This distribu¬ 
tion is very well described by an exponential distribu¬ 
tion at low redshift, and reasonably well in the range 
0.8 < z < 1.2, with an inverse slope —1/a increasing 
roughly as ^ (1 ± z)‘^-^. As a result, logic (/hds) 
roughly proportionally to — (1 ± independently of 

the threshold in H^a that is chosen (if H^a min> 4 A). 
However, we must add an important caveat to this result: 
the parameterization by Equations [21 and |2I is purely em¬ 
pirical, and constrained only by observations at ^ 0.1 

and 0.8 < 2: < 1.2. Any extrapolation to earlier epochs 
is hazardous. In particular, these equations suggest a 
never-ending increase of /hds with redshift, which is 
clearly not realistic when the universe is younger than 
a few 10^ years and galaxies cannot be in an HDS phase 
yet. 

5. MODELING AND INTERPRETATION 

To analyze the observed colors and spectral indices 
shown in the previo us figures, we use the Pegase.2 s pec- 
tral synthesis code (|Fioe fc Rocca-Volmeran^ll997D for 
broad-band colors and the 4000A break, together with 
its version at higher spectral resolut ion in the optical, 
Pegase-HR (|Le Borgne et all I2004D . needed to mea- 
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Fig. 9.— Distributions of the H(5 a index measured in four redshift 
bins. The slope measured on the SDSS sample (dotted line) for 
H(5 a^ 4 a is reported in the four panels. For GDDS galaxies, 
Monte-Carlo simulations were carried out from our sample of 25 
galaxies to estimate the real distribution at these redshifts; and the 
dashed line represents the best-fit pure expone ntia l distribution. 
The evolution of the slope a{z) is shown in Fig. UH top panel. 


sure Lick indices. The predictions of Lick indices with 
Pegase-HR and their accura cy with respect to othe r 
models are fully described in iLe Borgne et alJ 11200411 . 
The differences between the various models available to¬ 
day are small enough to reinforce the robustness of the 
results in this paper. 

In this section, we will compare the observed data with 
two classes of simple models; 
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Fig. 10.— Variations of /hds with the threshold H(5Amin for 
three redshift bins. This figure is directly derived from Fig El ex¬ 
cluding the bin 0.6 < 2 ; < 0.8 which does not contain enough 
galaxies to allow a reliable measurement of / hds - The evolution 
of the best-fit slope a'{z) is illustrated in Fig. 1111 bottom panel. 
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• “Quiescent” models, defined by smooth star for¬ 
mation histories, with continuous and modest star 
formation rates throughout the life of a galaxy. 

• “Burst” models, which correspond to galaxies that 
build up stellar mass in a few short and intense 
episodes of star formation. 

An archet ypal example of a quiescent model is t hat 
described by iLe Borgne &: Rocca-Volmerangd (|2002|) as 
doing a fairly good job of reproducing the optical and 
near-infrared colors of local Sc spiral galaxies. In this 
scenario a galaxy is fed by infall of zero-metallicity gas 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Redshift 

Fig. 11.— Evolutions of the exponential index a{z) of the 
H(5 a distribution (top panel) and of the exponential index a'{z) 
of /hds =f(H(5Amin) (bottom panel). The values of a and a' are 
reported from Fig.lSland llOl respectively. In the case of pure expo¬ 
nential distributions for H(5 a (equation the two panels should 
be identical. 


with a r = 8 Gyr e-folding time-scale. Inside the galaxy, 
the gas is simultaneously converted into stars following 
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a Schmidt law, with a specific efficiency = 0.1 Gyr“^ 
(1/z/ being the time required to convert one solar mass 
of gas into one solar mass of stars). A continuous star 
formation history can also be approximated by a series of 
small successive starbursts. The duty cycle of the bursts 
is encompassed by two evolutionary scenarios, hereafter 
referred to as “cycling”, or “recurring bursts”. In these 
models episodic events, which mimic constant star for¬ 
mation, occur during 100 Myr. These are followed by pe¬ 
riods of inactivity, lasting 100 Myr in one case (recurring 
bursts), and 0.9 Gyr in the other case (cycling). It can 
certainly be argued that our cycling model can hardly 
be considered to be continuous star formation, and in 
fact it might be better to think of the cycling model as 
representing a star-formation history that is intermedi¬ 
ate between the quiescent mode described here, and the 
“burst” mode described next. 

In contrast to our quiescent models, our burst models 
are intended to more closely reproduce the colors of lo¬ 
cal elliptical galaxies. In these models short time-scales 
(zy = r = 0.1 Gyr) are associated with the star-formation 
activity in the first three Gyr. After three Gyr, when 
the galaxy has converted most of the available gas into 
stars, galactic winds eject all the remaining gas from the 
galaxy, preventing further star formation. Two other, 
more extreme, examples of burst scenarios are a single 
instantaneous burst with solar metallicity, and a period 
of constant star formation truncated 0.5 Gyr after the 
collapse of the galaxy. 

Throughout the remainder of this paper we will set a 
redshift of z = 1.4 as the epoch of the last major star for¬ 
mation episode for all the burst-mode models considered. 
Glearly this redshift does not reflect the real redshift of 
formation for the galaxies being modeled, which we can¬ 
not strongly constrain on account of age-metallicity de¬ 
generacy. Instead 2 ; = 1.4 is simply a fiducial redshift 
that is supposed to mark roughly the end of the latest 
major event of active star formation that took place in 
the HDS massive galaxies. An instantaneous burst model 
with a redshift of formation z = 2 is also considered and 
shown in Figures 01 and to account for the old and 
red population of galaxies seen at ^ 1. In all these 

models, nebular emission lines are added to the stellar 
continuum. The fluxes in these lines are computed with 
Pegase.2 for a mean H II region consistently with the 
on-going and recent star-formation activity in the galaxy. 
When the star-formation activity is intense, these lines 
fill-in the absorption stellar lines and often lead to very 
negative values for the [Oil] and US equivalent widths. 

Tracks corresponding to a representative set of quies¬ 
cent and burst models are superposed on Figures IHHZI 
The tracks shown on Figure 0 are particularly illumi¬ 
nating, as these suggest an evolutionary link between 
the z = 1 massive HDS galaxies and the local massive 
and red galaxies. Indeed, the HDS galaxies that are al¬ 
ready massive at 2 : = 1 can only be identified as the 
progenitors of some local massive galaxies: they cannot 
lose much of their stellar mass (which can be caused by 
the death of a fraction of the stellar population) during 
their subsequent evolution. Therefore, the distant HDS 
galaxies necessarily must be linked to the present mas¬ 
sive galaxies, most of which are passively evolving. It is 
clear that this evolution can only he modeled by a break in 
the star-formation history: the z 1 HDS galaxies are 


just leaving a burst mode with intense star formation 
activity (suggested by their already large stellar masses) 
which ended 0.5 — 1 Gyr before the epoch of observa¬ 
tion, at ^ ^ 1.2 —1.5. When seen 1-2 Gyr after the burst, 
at 2 : 0.8, they are evolving almost passively with little 

star formation, as described in more details below. 

We emphasize that the bulk of the massive SDSS galax¬ 
ies seen at redshift 2 : = 0.05 — 0.1 have spectral fea¬ 
tures characteristic of evolved galaxies with very little 
star formation activity. Their small H(5 a (Fig. |n|), large 
D^4000 and red g — r color (Fig. ^ are all well mod¬ 
eled by a passively evolving stellar population, as il¬ 
lustrated by the models shown in Figures IHHZI High- 
redshift analogs to these objects are seen at ^ ^ 1 in 
the population of GDDS galaxies marked by solid cir¬ 
cles in figures 01 These objects present the same small 
US A values, large D^4000 indices, and red g — r colors 
as evolved galaxies in the local universe^^. These objects 
connect directly onto a number of recent studies of mas¬ 
sive passi vely evolving galaxies out to even higher red- 
shifts ('e.g.lMcCarthv et all2004lGlazebrook et all2004 


A comparison of the “quiescent” models (which all 
have a smoother star-formation history than the burst 
models) with the data is also interesting. Figure 01 shows 
that continuous star-formation for a long period of time 
(Sc model) leads to intermediate values of the RSa index 


(2 - 3 A). This is not true if the star-formation is in¬ 
termittent: when the star-formation episodes are spaced 
by more than 100 Myr (which is the case of our recur¬ 
ring and cycle models), the H^a index shows substantial 
variations from very positive (if star-formation is halted 
abruptly) to very negative values (because of nebular 
emission lines). The cycle and recurring burst models 
shown in Figure 0 also make it clear that if the star for¬ 
mation episodes are too frequent (i.e. repeated on pe¬ 
riods smaller than a very few Gyrs), the large D^dOOO 
indices and red rest-frame colors of the massive galaxies 
at ^ = 0 or 2 : = 1 can never be reached. 

A particularly nice aspect of modeling RSa and 
Dn4000 is that models used for interpretation of these 
two quantities are sufficiently robust that uncertainties 
in them are unlikely to lead to wild misinterpretations 
of the data. Indeed, we have checked that a change in 
the IMF, within reasonable limits, does not affect the 
predicted USa index by more than 2 A and D^4000 by 
more than 0.2. To some extent, the uncertainties in the 
stellar evolution tracks can be explored by looking at the 
differences between the predi ction of the Pegase mod ¬ 
els and the Galaxev m odels (|Brnznal Gharl^l2003ll . 
iLe Boreme et al.l (|2004ll showed that the differences in 
the predicted H^a indices are smaller than 1 A at any 
metallicity. We also verified that the differences on the 
predicted D^dOGO index are smaller than the errors in 
the data. However, uncertainties on the effect of the 
non-solar abundances of a elements remain unclear and 
might be responsible for slight deviations of the model 
tracks from the data. The general trends will, however, 
remain robust. 


Note that the stellar mass cutoff we chose when defining our 
sample corresponds (roughly) to the pivot mass in the observed 
bimodality between blue and red galaxies in the local universe 
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It is worth emphasizing that strong 116 in absorption 
is, strictly speaking, a signature of a recent break in star- 
formation, rather than of any preceding star-burst ac¬ 
tivity. Moreover, a complete halt in the star-formation 
activity is not needed to obtain an HDS system. It is 
the contrast between the current and recent-past star- 
formation rates that must be large. 

To illustrate the importance of this contrast, we can 
show that the current star-formation rate is small in com¬ 
parison to its average past value by comparing the char¬ 
acteristic growth time-scale of a galaxy to the age of the 
universe at its redshift. Following Paper V, we define 
this growth time-scale (expressed in Gyr) as the ratio 
M^/SFR\ it is the time needed for a galaxy to assemble 
its current mass if its SFR is fixed at its current value. If 
the growth time-scale of a galaxy is larger than the age 
of the universe at its redshift, we can conclude that such 
a galaxy is in a mode of passive evolution (relative to its 
past star-formation history). The SFR values of our high 
redshift galaxies, estimated from the [OIIJA3727 emission 
line and assuming Ay=I, are listed in Tabled The typ¬ 
ical values of the on-going SFR are 1-2 Mq yr“^ for HDS 
galaxies (and O.I — I Mq yr“^ for the non-HDS, “older” 
galaxies). Although these values are non negligible, they 
are still small with respect to the past mean SFR: Tabled 
shows that M^/SFR is larger than the age of the uni¬ 
verse, in a ratio of 3 for the HDS galaxies, and 10 for 
the non-HDS galaxies. This is in good agreement with 
the results of Paper V, and it means that the massive 
galaxies are already in a mode of passive evolution at 
z = 1. It also illustrates the large contrast between the 
recent (:^ I Gyr ago) large bursts of star formation and 
the modest values of the on-going SFR in HDS systems. 

The classical way to study “post-starburst” galaxies 
(i.e. with completely truncated star formation) is to 
focus on the sub-sample of HDS galaxies with little 
(or no) emission lines, as measured by their equivalent 
widths. The most well-known categorization of these 

S galaxies is the ‘E+A’ population llpressler fc GunrJIl^ 

h et alJlQQalPoggianti et alJl999l:lZabludofF et alJ 

, with strong RS absorption lines coupled with lit¬ 
tle (or no) [OH] emission. To facilitate comparison with 
previous work, we conclude this section by investigating 
the fraction of ‘E+A’ galaxies in our HDS population. 
Table Q includes measurements of [OH] equivalent width, 
with emission lines represented by negative values^^. As 
described below, our GDDS sample contains one or two 
E+A galaxies, depending on how error bars are treated. 

Figure IT^ shows the distribution of the two samples of 
massive galaxies on a EW(H(5)-EW[0 II] diagram. The 
SDSS sample is represented with the same color-coding of 
the number density as in previous figures. The individual 
GDDS galaxies are shown with their error bars. Some 
very simple closed-box evolutionary models computed 
with Pegase are over-plotted: these models have expo¬ 
nentially decreasing star-formation rates with e-folding 
time-scales r = 3,10,100, 200 and 1000 Myr (Column I 

Note that, for this analysis, we chose to use the definition 
of the equivalent width EW(H(5) given in Table |2| instead of the 
Lick H(5 a index used earlier in this paper in order to make our re¬ 
sults directly comparable to previou s work ([^braha m et iDIml 
IBalogh et al.ll99ai(;oto et all200aiTran et al.l2004l). When mea¬ 

sured on the SDSS sample, the average value of EW(H(5)—H(5 a is 
0.50 A, with an rms dispersion of 0.78 A. 



Fig. 12. — Selection of E+A galaxies (defined by the dashed box: 
EW(H(5)> 4 A and EW[0 II] > —5 A). The shaded colored region 
corresponds to the SDSS sample (same color-coding as in Fig. 1!^. 
The individual GDDS galaxies of our sample are shown on the top 
of the SDSS sample, with their error bars. The evolutions predicted 
by closed-box exponential models (solid lines, with arrows marking 
every doubling in age) are given for various e-folding time-scales r. 
The fraction of galaxies in the E+A phase (or simply in the HDS 
phase with EW(H(5)> 4 A) is much larger in the GDDS sample 
than at low redshift. Moreover, only models corresponding to a 
rapid decrease of the star-formation rate (r < 300 Myr) reach the 
E+A phase (see also Table EJ- 


TABLE 3 

Closed-box models with exponentially decreasing SFR. 


byr) 

EW(H5) 

(A) 

Age 

(Gyr) 

AtHDS AtE+A 
(Gyr) (Gyr) 

0.003 

10.0 

0.25 

1.56 

1.56 

0.01 

10.0 

0.30 

1.54 

1.53 

0.1 

8.8 

0.60 

1.48 

1.11 

0.2 

7.1 

0.80 

1.51 

0.67 

1.0 

3.5 

1.80 

0.00 

0.00 


of Table IS)). All these models form stars at solar metal- 
licity. The arrows give the direction of the evolution 
in the diagram: EW[0 II] is highly negative soon af¬ 
ter the maximum of the star formation because of the 
presence of strong emission lines, and increases contin¬ 
uously afterward. The models with r < 300 Myr pass 
through an HDS phase (if defined by EW(HJ)> 4 A) 
and an E+A phase (if defined by an HDS galaxy with 
EW[0 II] > —5 A) which is illustrated by the dashed line 
in Figure IT^ 

The maximum EW(HJ) values reached by our mod¬ 
els, and the corresponding ages, are given in Table 01 
Golumns 2 and 3. This table also records (Golumns 
4 and 5) the total time spent in the HDS and E+A 
phases. The table shows that (i) the maximum EW(H(5) 
is reached sooner for rapidly decreasing SERs, (ii) the 
time spent in the HDS phase is almost independent of 
the e-folding time-scale when this phase is reached, and 
(iii) a given galaxy spends less time in the E+A phase 
when the SER declines more slowly. On this basis, it is 
clear from Eigure^lthat the proportion of galaxies show- 
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ing signs of recently truncated star formation (r < 200 
Myr and EW(H^)> 4 A) is higher for the GDDS sample 
than for the SDSS sample. 

Models in Table |3| and Figure IT^ can be used to make 
a theoretical prediction regarding the evolution of E+A 
systems: while models with r < 300 Myr have both an 
HDS and an E+A phase, all models that have an HDS 
phase do not necessarily have a E+A phase. (For in¬ 
stance, models with 0.3 < r < 1.0 Gyr). Therefore, 
depending on how fast the star-formation is halted, one 
might expect that the evolution of the fraction of E+A 
galaxies would show milder evolution than the fraction 
of HDS galaxies. Strictly speaking, 4^3% (1/25) of the 
GDDS galaxies are E+A at ^ = 0.6 — 1.2 and 3.1% 
(1561/50,255) of the SDSS sample at ^ = 0.05 — 0.1 meet 
the E+A definition. (The GDDS galaxy 02-1543 might 
be counted as an additional E+A, depending on how one 
treats its error bar.) The fractions at low and high red- 
shifts being roughly compatible, we don’t observe any 
evolution in the fraction of E+A galaxies. However, the 
small number statistics involved and our large error bars 
leave room for further studies on this topic. 

It is of interest to note that a number of recent pa¬ 
pers derive a (1 + 2 :)"^ evolution to describe the rise in 
the merger rate as a function of redshift, with m indices 
that are similar to the exponential index we obtain for 
the rise in the fraction of HDS systems as a function 
of redshift: compare our va lue of m = 2.5 + 0.7 with 
m = 3.4 ± 0.6 obtained bv iLe Fevre et and 

m = 2.5±0.7 obtained bv iPatton et al.l (|2002ll . However, 
the evolution in the merger rate is poorly constrained. 
For example, l arger values ye al so obtained, such as 
m = 4 ± 1 bv iRnshetnikoyl (l2000fl and m = 4 — 6 ac¬ 
cording to lGonselice et al.l(|2003D . At this stage we think 
the most that can be said is that our results may hint 
at a relationship between merging and the HDS phase. 
A similar re lationship is suggested by the recent work of 
iGotol (|2005ll on 266 E+A galaxies in the SDSS. If a link 
exists between HDS galaxies and mergers, the morpholo¬ 
gies of HDS galaxies should reflect these dramatic events. 
Figure m shows that many of the HDS galaxies (left col¬ 
umn) indeed have disturbed morphologies which can be 
signs of merging. Moreover, the ratio between irregular 
and relaxed morphologies should roughly scale as the ra¬ 
tio of time-scales for merging (:^ 0.5 Gyr) and spectral 
evolution 0.3 — 2 Gyr). This statement seems gener¬ 
ally consistent with Figured which shows several HDS 
galaxies with almost relaxed morphologies (02-1724, 02- 
1777, 22-0315), but further discussion of this topic is re¬ 
served for an upcoming paper in this series (Abraham et 
al. 2006, in preparation). 

6. CONCLUSIONS 

By comparing two mass-selected samples of galaxies at 
low (SDSS) and high (GDDS) redshift, we find that the 
fraction of galaxies with strong HJ absorption lines in¬ 
creases with redshift up to 2 : = 1.2. This evolution can 
be parameterized as a function of redshift: (1 + 

It does not depend on the threshold chosen for 116a to 
define HDS galaxies, and it reflects a similar evolution for 
the distributions of the 116 equivalent widths. By model¬ 
ing this variation using Pegase models as a function of 
D^4000 and the rest-frame g — r color, we conclude that 
this evolution originates from a fundamental change in 


the mode of star formation in massive galaxies that oc¬ 
curred somewhere around 2 : 1.4 — 2. Our models sug¬ 

gest that the star formation rates drop sharply in most 
massive galaxies at that epoch, which is consistent with 
the measures of star formation rate densities presented 
in Paper V. 

This drop in the star formation rate, together with 
the already large stellar masses of these 2 : 1 galax¬ 

ies, suggests very active phases in the star-formation at 
higher redshifts. The massive obscured starbursts being 
detected with the Spitzer Space Telescope at ^ 2 — 2.5 

could very well be the progenitors of our sample of 2 : 1 

galaxies: the “old” galaxies seen in our sample at 2 : = 1, 
with small H^a indices and large D^4000 breaks, have 
necessarily experienced their last significant episodes of 
star-formation at 2 : > 1.8 and assembled most of their 
stellar mass before that time. These episodes must have 
been very intense (at least 30 MQyr~^ on average), 
given the large stellar mass of these galaxies and the age 
of the universe at this redshift. Moreover, some of the 
^ = 2 — 2.5 starbursts could also have spanned larger 
timescales (2 — 3 Gyr) in which case they could become 
HDS galaxies at 2 : ^ 1, on the condition that their star- 
formation declined at ^ ^ 1.4, on time-scales shorter 
than a few hundred Myrs. This pictu re is in fairly good 
agreement with the recent findings oflTreu et all II2005D 
who estimate that the massive field spheroidal galaxies 
seen at 2 : = 0.2 — 1.2 in the GOODS-N field formed most 
of their stellar mass at 2 : > 2, with subsequent activity 
continuing to lower redshifts. 

Finally, our study of spectral features in massive galax¬ 
ies from ^ = 1.2 to the present shows that the redshift 
range 2 : = 1.4 — 2 corresponds to the epoch at which the 
star-formation in these galaxies has changed from an ac¬ 
tive phase to a passive one. An echo of this change in 
the mode of star-formation activity is imprinted in the 
strong H(5 absorption lines seen in massive galaxies at 
z ^ 1. 

7. ACKNOWLEDGMENTS 

The authors wish to thank Jarle Brinchmann for many 
stimulating discussions and for sharing with us a num¬ 
ber of plots which motivated much of this work. We are 
also in debt to the anonymous referee whose comments 
helped to improve the quality of this paper. Observa¬ 
tions were obtained at the Gemini Observatory, which is 
operated by AURA Inc., under a cooperative agreement 
with the NSF on behalf of the Gemini partnership: the 
NSF (US), PPARG (UK), NRG (Ganada), GONIGYT 
(Ghile), ARG (Australia), GNPq (Brazil) and GONIGET 
(Argentina) and at the Las Gampanas Observatory of 
the OGIW. RGA acknowledges support from NSERG 
and the Province of Ontario. KG and SS acknowledge 
support from the David and Lucille Packard Founda¬ 
tion. Funding for the creation and distribution of the 
SDSS Archive has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National 
Aeronautics and Space Administration, the National Sci¬ 
ence Foundation, the U.S. Department of Energy, the 
Japanese Monbukagakusho, and the Max Planck Soci¬ 
ety. The SDSS Web site is http://www.sdss.org/ The 
SDSS is managed by the Astrophysical Research Gon- 
sortium (ARG) for the Participating Institutions. The 
Participating Institutions are The University of Ghicago, 












16 


D. Le Borgne et al. 


Fermilab, the Institute for Advanced Study, the Japan 
Participation Group, The Johns Hopkins University, the 
Korean Scientist Group, Los Alamos National Labora¬ 
tory, the Max-Planck-Institute for Astronomy (MPIA), 
the Max-Planck-Institute for Astrophysics (MPA), New 


Mexico State University, University of Pittsburgh, Uni¬ 
versity of Portsmouth, Princeton University, the United 
States Naval Observatory, and the University of Wash¬ 
ington. 


REFERENCES 


Abazajian, K., Adelman-McCarthy, J. K., Agiieros, M. A., et al. 
2004, AJ, 128, 502 

Abraham, R. G., Glazebrook, K., McCarthy, P. J., et al. 2004, AJ, 
127, 2455 [Paper I] 

Abraham, R. G., Smecker-Hane, T. A., Hutchings, J. B., et al. 

1996, ApJ, 471, 694 

Abraham, R. G. et al. 2006, AJ, in prep. 

Baldry, I. K. & Glazebrook, K. 2003, ApJ, 593, 258 
Baldry, I. K., Glazebrook, K., Brinkmann, J., et al. 2004, ApJ, 600, 
681 

Balogh, M. L., Morris, S. L., Yee, H. K. C., Carlberg, R. G., & 
Ellingson, E. 1999, ApJ, 527, 54 
Bell, E. F., et al. 2004, ApJ, 608, 752 
Blanton, M. R., et al. 2003, ApJ, 592, 819 

Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, MNRAS, 
351, 1151 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 
Cardiel, N., Gorgas, J., Cenarro, J., & Gonzalez, J. J. 1998, A&AS, 
127, 597 

Chen, H., McCarthy, P. J., Marzke, R. O., et al. 2002, ApJ, 570, 
54 

Cimatti, A., Daddi, E., Renzini, A., et al. 2004, Nature, 430, 184 
Cole, S., et al. 2001, MNRAS, 326, 255 

Conselice, C. J., Bershady, M. A., Dickinson, M., & Papovich, C. 
2003, AJ, 126, 1183 

Couch, W. J. & Sharpies, R. M. 1987, MNRAS, 229, 423 
Dressier, A. & Gunn, J. E. 1983, ApJ, 270, 7 
Dressier, A. & Shectman, S. A. 1987, AJ, 94, 899 
Fasano, G., &: Franceschini, A. 1987, MNRAS, 225, 155 
Fioc, M. & Rocca-Volmerange, B. 1997, A&A, 326, 950 
—. 1999, astro-ph/9912179 

Firth, A. E., Somerville, R. S., McMahon, R. G., et al. 2002, 
MNRAS, 332, 617 

Fontana, A., Pozzetti, L., Donnarumma, I., et al. 2004, A&A, 424, 
23 

Gehrels, N. 1986, ApJ, 303, 336 

Glazebrook, K., Abraham, R. G., McCarthy, P. J., et al. 2004, 
Nature, 430, 181 [Paper III] 

Goto, T. 2005, MNRAS, 357, 937 

Goto, T., Nichol, R. C., Okamura, S., et al. 2003, PASJ, 55, 771 
Juneau, S., Glazebrook, K., Grampton, D., et al. 2005, ApJ, 619, 
L135 [Paper V] 


Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, 
MNRAS, 341, 33 

Kewley, L. J., Geller, M. J., & Jansen, R. A. 2004, AJ, 127, 2002 
Kroupa, P. 2001, MNRAS, 322, 231 

Le Borgne, D. & Rocca-Volmerange, B. 2002, A&A, 386, 446 
Le Borgne, D., Rocca-Volmerange, B., Prugniel, P., et al. 2004, 
A&A, 425, 881 

Le Fevre, O., Abraham, R., Lilly, S. J., et al. 2000, MNRAS, 311, 
565 

Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 
McCarthy, P. J., Carlberg, R. G., Chen, H.-W., et al. 2001, ApJ, 
560, L131 

McCarthy, P. J., Le Borgne, D., Crampton, D., et al. 2004, ApJ, 
614, L9 [Paper IV] 

Panter, B., Heavens, A. F., & Jimenez, R. 2004, MNRAS, 458 
Paterno, M. 2004, Calculating Efficiencies and 
their Uncertainties, Tech. Rep. FERMILAB-TM- 
2286-CD, Fermilab Technical Memo, see also 
http: //home. f nal. gov/^paterno/images/ef f ic . pdf 
Patton, D. R., Pritchet, C. J., Carlberg, R. G., et al. 2002, ApJ, 
565, 208 

Poggianti, B. M., Small, L, Dressier, A., et al. 1999, ApJ, 518, 576 
Quintero, A. D., Hogg, D. W., Blanton, M. R., et al. 2004, ApJ, 
602, 190 

Reshetnikov, V. P. 2000, A&A, 353, 92 
Salpeter, E. E. 1955, ApJ, 121, 161 
Schmidt, M. 1968, ApJ, 151, 393 

Shioya, Y., Bekki, K., & Couch, W. J. 2004, ApJ, 601, 654 
Tran, K. H., Franx, M., Illingworth, G. D., et al. 2004, ApJ, 609, 
683 

Treu, T., et al. 2005, ApJ, 633, 174 
Worthey, G. & Ottaviani, D. L. 1997, ApJS, 111, 377 
Yang, Y., Zabludoff, A. I., Zaritsky, D., Lauer, T. R., & Mihos, 
J. C. 2004, ApJ, 607, 258 

Zabludoff, A. 1. 1999, in lAU Symp. 186: Galaxy Interactions at 
Low and High Redshift, 125 

Zabludoff, A. L, Zaritsky, D., Lin, H., et al. 1996, ApJ, 466, 104 


